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Abstract. We present a high spatial and spectral resolu- 
tion 2-D echelle spectrogram of the Narrow-Line Region 
in the Seyfert 2 galaxy NGC 1386. This Seyfert galaxy was 
observed with CASPEC in the wavelength range 5270- 
7725 A which covers the Ha and the [Nil] lines. With 
the use of spatially high resolved images taken with the 
WFPC2 aboard the Hubble Space Telescope we could 
identify individual components of the Narrow-Line Region 
in our spectra. A Gaussian decomposition of the spectra 
revealed 9 distinct emission-line complexes. The bright- 
est component is blue-shifted by —120 ± lOkms"^ with 
respect to the systemic velocity and shows an offset of 
« — 1'.'6 relative to the nucleus of the galaxy. The true 
nucleus of NGC 1386 has a much lower apparent Ha lu- 
minosity than this component. The nucleus is probably 
highly absorbed. Although the majority of the Narrow- 
Line Region components follows a regular velocity field, 
we find evidence for a separate kinematic component. The 
Narrow-Line Region is aligned anti-parallel to the radio- 
jet which propagates from the center of NGC 1386 to the 
south. 

Key words: galaxies: active - galaxies: Seyfert - galaxies: 
individual: NGC 1386 - galaxies: kinematics and dynamics 
- galaxies: structure 



velocity components, such as inflow or outflow motions, 
may be superimposed onto the regular rotational velocity 
field (e.g., Osterbrock 1993). The morphology of the ion- 
ized gas component often coincides with radio structures, 
seen in VLA maps (e.g., Ulvestad & Wilson 1984; Pogge 



1989). Such extended radio structures (e.g., radio jets) are 
often aligned with the optical emission (cf. Capetti et al. 
Gallimore et al. 



1997 



1996) 



First hints for a complex structure of the emission line 
region of A GN w ere mentioned by Walker (quoted by Bur- 
bidge et al. [1959D in NGC 4151 und NGC 1068. In the late 
60s and later on particularly for NGC 1068 strong indi- 
cation for a structured emission line region und indiviual 
substructures were found (Walker 1968|; An derson 1970 



Glaspey et al. 1976 a, b; AUoin et al. 1983| ). Recent in- 



vestigations of nearby Seyfert galaxies with high spatial 
resolution, as it is achieved with the Hubble Space Tele- 
scope (HST), showed that it is feasible to resolve indi- 
vidual NLR components in the nearest Seyfert galaxies 
(e.g., NGC 1068, Macchetto et al. |1994| ; NGC 4151, Winge 
et al. |1997t Mkn3, Mkn78, Mkn348, Capetti et al 



1995 



1996b) on scales as small as O'.'l. In the case of NGC 1068 



the morphology of the NLR can be characterized as a 
clumpy structure of several dozen individual emission line 
knots. 



1. Introduction 

The emission line regions of active galactic nuclei (AGN) 
can be divided into the Broad-Line Region (BLR) and the 
Narrow-Line Region (NLR). The study of velocity fields 
of the narrow emission-line region in nearby AGN offers 
an exciting perspective to study details of the circumnu- 
clear kinematics. Besides the primary component, which 
is given by the rotation of the galaxy, additional velocity 
components are likely to be present. Especially, in AGN 

Send offprint requests to: jrossa@astro.ruhr-uni-bochum.de 

* Based on observations collected at the European Southern 
Observatory, La Silla, Chile 
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1.1. NGC 1386 

The galaxy NGC 1386 (Fig. [l]) has been cla ssified as a 
Seyfert type 2 galaxy by Phillips & Frogel ( |l980| ). The 



spectra show characteristic strong emission lines of Ha, 
H/3, [Nil], and [Oiii] with expected line ratios and line 
widths typical for a Seyfert nucleus. 

There is some debate on the morphological galaxy 
type, whether it is Sa or SO, due to the relatively high 
inclination of i « 71°. TuUy ( |1988D hsts it as SO while 
generally a Sa classification is assumed (e.g., Sandage & 
Tammann 1981; Weaver et al. |1991| ). The assumption from 
the Weaver et al. study was based on the line ratios found 
for the outer regions of NGC 1386 which were characteris- 
tic of H II regions. In the recent literature one can also find 
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Table 1. Basic data on NGC 1386 



Galaxy 


R. A. (J 2000)" 


Dec. (J 2000) 


Hubble typef 


AGN type 


D WHi"^ a X b i mji^ 


NGC 1386 


03''36'"46=4 


-35°59'58'.'l 


Sa 


Sy2 


12.2Mpc 918kms-i 3.'5 x l!3 71° 11.24 



" All data have been taken or have been calculated from the RC3 (de Vaucouleurs et al. 1991), except where indicated 
' taken from Weaver et al. (1991) 



taken from Huchtmeier & Richter (1989) 

taken from the NASA Extragalactic Database (NED) 



classifications with SBO+ (Tsvetanov & Petrosian 1995| ), 
although there are no clear indications for the existence 
of a bar. 




Fig. 1. The Seyfert 2 galaxy NGC1386. This is the cen- 
tral region of the HST archive image obtained with the 
WFPC 2, and was taken with the F606W filter. North is 
up, East is to the left. The dust lanes, which run across 
NGC 1386, are easily recognizable. 



Since attention was first paid on this galaxy, there was 
some uncertainty whether NGC 1386 belongs to th e For - 
nax galaxy cluster or not. Huchtmeier & Richter ( 1989 ) 
list a heliocentric corrected systemic velocity of 918 ± 
34kms^^, derived from Hi measurements which is in 
agreement wit h the value of 890 ± 20 kms~^ obtained by 
Weaver et al. (1991), as derived from the rotation curve. 

We assume a Hubble parameter, that will be used 
throughout this paper, of Hq = 75kms~^ Mpc~^. Assum- 
ing no effect from the Fornax Cluster the radial veloc- 



ity, obtained by Huchtmeier & Richter ( |1989|) , translates 
to a distance to NGC 1386 of 12.2 Mpc. Basic data on 
NGC 1386 are listed in Table 0. 

The distance to the Fornax Cluster has been deter- 
mined by McMillan et al. (1993), using the Planetary 
Nebula Luminosity Function (PNLF). From their study of 
three member galaxies of the Fornax Cluster they derived 
a distance of 16.9 ± 1.1 Mpc. This value is coincident with 



the value that TuUy (1988) listed in his Nearby Galaxies 
Catalog. More recent measurements using Cepheid based 
distances yielded 18.6 ± 0.6 Mpc (Madore et al. WdM). 



They also derived a mean velocity of the Fornax Cluster 
oi V = 1441 km s~^ with a = ±342kms~^. We therefore 
assume that NGC 1386 is located in the foreground of the 
Fornax Cluster. 



Being one of the closest Seyfert nuclei, NGC 1386 is 
an ideal candidate for high spectral and spatial resolu- 
tion investigations of the NLR, and we therefore chose 
NGC 1386 for a detailed investigation of the NLR kine- 
matics and morphology. One of the techniques applied to 
this object class besides narrowband imaging (e.g., Evans 
et al. 119911; Lynds et al. 119911; Pogge & De Robertis f993|) 



is the speckle masking method (e.g.,Mauder et al. 1992) 
although there are alternative techniques that have al 
ready been applied to the nearest Seyfer t galaxies success- 
full y (e.g ., Wagner & Appenzeller 1988; Dietrich & Wag- 



ner [19981) . However, the first attempt to resolve individ- 
ual NLR components within NGC 1386 was mad e using 
speckle masking observations (Mauder et al. 1992 ). From 
their study it was inferred that the morphology of the gas 
distribution is cone-like. 



Weaver et al. ( |199lD studied the NLR in NGC 1386 
with medium spectral resolution. They reported on the 
velocity field as a combination of a normally rotating com- 
ponent and a component undergoing high velocity infall 
or outflow, as derived from the shape and width of the line 
profiles. The line width amounts to 540 ±40 km s~^ at the 
nucleus with a double peak line profile. These measure- 
ments are in agreement with our better resolved line pro- 
files (due to higher spectral and spatial resolution) which 
yield a FWHM of 470 ± 20km 8"^ Furthermore, they re- 
ported that the optical nucleus is displaced from the kine- 
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matical center of the rotation curve by 1" — 1'.'5 which 
corresponds to 60-90 pc. 

2. Observations and data reduction 

2.1. Echelle spectroscopy 

The Seyfert galaxy NGC 1386 was observed at ESQ La 
Silla, Chile (Nov. 16, 1992) with the Cassegrain Echelle 
Spectrograph (CASPEC) attached to the ESQ 3.6m tele- 
scope. The achieved spectral resolution with the short 
camera is R = ^ « 18000 which corresponds to Aw = 
17kms~^. With the use of the red cross-disperser the 
echelle spectrum covers the range 5270-7725 A. The slit 
covered an area of 8" x 1" and was oriented at p. a. 16°. 
The CCD chip, a TEK512M-12 with 512 x 512 pixels 
of 27/impix~^ offers a spatial scale of 0'.'65pix^^. The 
NGC 1386 spectrum had an integration time of 5400 sec. 

The data reduction was performed with the IRAF[| and 
MIDASQ software packages. The spectrum was bias cor- 
rected creating a master bias which was subtracted from 
the spectrum. The tracing was determined from standard 
star exposures. In this process a Chebychev polynomial of 
3'''^ order was fitted along the dispersion axis of each aper- 
ture. The master trace was defined by the spectrum with 
the smallest rms deviations. The echelle spectrum was 
then background corrected making use of the interorder 
sections of the 2D echelle spectrogram. 

The fiatfield correction was done in two different steps. 
Due to spatial offsets and different slit lengths, we had 
to create several groups of flatfields. The internal scat- 
ter of the spatial offsets was less than 1 pix within a flat- 
field group. The deviation between the individual flatfield 
groups had offsets of the order of Ay w 2 pix. For each 
group a normalized flatfield was computed which were 
used to correct for overall sensitivity variations of the cor- 
responding science frames. Due to the flexure of the spec- 
trograph the echelle object spectra had small offsets with 
respect to the corresponding flatfield exposure. Hence, we 
were forced to restrict the 2D fiatfield correction to the 
inner 3"while the 1'.'5 wide stripes above and below this 
region were not corrected in this step. The normalized 
fiatfield was set to unity for the region of these stripes. 
The criterium for the selection of these scans based on an 
intensity of less than w50 % of the maximal intensity on 
the position of the individual scan. In the second step the 
final extracted ID spectra of the inner 6"of the individual 
spatial scans were fiatfield corrected with a ID fiatfield 
generated. It was created by normalizing the individual 
scans of the used fiatfield for the 2D correction with the 
average ID flatfield of it. Finally we are able to make use 



^ IRAF is distributed by the National Optical Astronomy 
Observatories, which is operated by the Association of Univer- 
sities for Research in Astronomy, Inc. (AURA) under cooper- 
ative agreement with the National Science Foundation. 

^ Munich Image Data Analysis System, distributed by ESO 



of the 2D echelle spectrograms covering the NLR gas of 
NGC 1386 out to a radius of « 3" from the nucleus. 

Subsequently the data were wavelength calibrated by 
using a ThAr lamp reference spectrum. The unit step was 
set to 0.1464 A. The scatter of our wavelength calibration 
had an rms of 0.014 A. The spectral resolution of our spec- 
tra was obtained from measurements of several night-sky 
lines in our spectra and was found to be 0.3 A. A final 
flux calibration was applied using the standard HR3454 
(T^Hya) which was reduced in the same way as the 2-D 
spectrum of NGC 1386. The extracted spectrum was fi- 
nally merged and a two-dimensional image was created 
from the 11 pixel scans. To correct for the effect of slit illu- 
mination, we divided our 2-D galaxy echelle-spectrogram 
by an identically processed normalized sky-fiatfield ex- 
posure. Finally a cosmics correction was applied to the 
spectrum using the MIDAS modily/gcurs routine. 

In Fig. H we present the spectrum of NGC 1386 from 
the central 2'.'6. The typical strong emission lines of Seyfert 
galaxies, e.g.. Ha, [Nil], [Sii] are visible. Above 6900 A the 
spectrum is dominated by strong atmospheric absorption 
bands (b-band, A-band). 



1 1 1 1 1 1 1 1 1 1 1 1 

NGC13B6 
p. a, 16° 









5500 



6000 6500 
Wavelength \k] 



7000 



Fig. 2. The averaged, central 2'.'6 of the NGC 1386 2D 
echelle spectrogram are shown. The fiux in units of 
10^^^ erg s^^ cm^^ A^^ is plotted versus the wavelength 
in the range of 5270-7200 A. The interstellar NaD ab- 
sorption line at A « 5900 A is clearly visible. 



2.2. Imaging 

We retrieved high resolution images of NGC 1386 from the 
Space Telescope Science Institute (STScI) data archive 
(from the survey by Malkan et al. 1998, Prop. -No. 5479). 



This image had been obtained with the WFPC 2 aboard 
the Hubble Space Telescope (HST) . It was taken with the 
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F606W filter (A„ean=5843.3 A, FWH M^15 78.7 A). Fur- 
ther details are given in Biretta et al. (1996). 

Additionally, further WFPC 2 images of NGC 1386 be- 
came available from a survey of Seyfert galaxies carried 
out by Wilson et al. (Prop. -No. 6419). These images had 
been taken in two narrowband and two broadband filters. 
In our study we made use of the F547M and F658N fil- 
ter images. The narrowband (F658N) image covers the 
HQ!-h[Nn] emission within the filter FWHM of 28.5 A. The 
WFPC 2 images of NGC 1386 were used to identify the 
individual components of the NLR in our 2-D spectrum 
which will be described in detail in Sect. 4.3. 



with spatial information as a function of velocity. In Fig. |^ 
the spatial scans of the velocity space transformed spec- 
trum in the [N ll] A6583 line are shown. In order to make 
a quantitative approach, a decomposition using elliptical 
Gaussians was applied to the spectrum. This yielded the 
position and velocity of the individual components in the 
I — u-map as well as the line widths for each individ- 
ual cloud in the / — v space. The detailed description of 
the process of decomposing a 2-D echelle spectrogram is 
given by Dietrich & Wagner ( 1998 ) for the Seyfert galaxy 
NGC 1068. 



3. Analysis 



In addition to the pure emission line fluxes the calibrated 
spectrum also contains the radiation from the underly- 
ing stellar population (continuum). To remove this ef- 
fect in a simple approach we calculated a linear contin- 
uum fit to our spectrum across the emission lines of in- 
terest (e.g., Ha, [Nil]). This linear interpolation which 
based on the mean flux in continuum windows («10A 
width) at either side of the Hq:+ [N ii] doublet emission line 
complex (Ag;„^=6538A, A5;^rf=6675A) was determined for 
each spatial scan and subtracted. 

The 2-D echelle spectrogram was transformed from 
wavelength space into velocity space as described by Di- 
etrich & Wagner ( 199^ ) . We used the heliocentric velocity 
v^ys = 918 ± 34kni£^of NGC 1386, as given by Hucht- 



meier & Richtcr ( 1989 ). We checked this value by mea- 
suring the centroid of the redshifted interstellar Nal line 
and applying the appropriate heliocentric correction. This 
results in Vsys — 924 ± 51 kms~^, in good agreement with 



Huchtmeier & Richter ( |1989D . 

Since our spectrum covers the wavelength range 5270- 
7725A, we have several strong emission lines which are 
useful for the kineniatical study of the line emitting gas. 
We use Ha and the [Nil] lines for our kinematical in- 
vestigation. The less intense [S ii] AA 6717,6731 doublet, 
[Ol] A6300 and the [Fevii] AA 5721,6087 lines wiU be used 
for the investigation of the spatial variations of the physi- 
cal conditions in the individual NLR emission-line clumps. 




-300 -200 -100 
V [km/s] 



-300 -200 -100 100 
V [km/s] 



Fig. 3. The spatial profiles of the [Nil] A 6583 A line are 
shown. The velocity relative to the systemic velocity is 
plotted versus the flux in units of 10~^^ ergs~^ cm~^ A 
The left panel shows 6 southern scans plus the center 
scan, each separated by 0'.'5 from the nucleus, and the 
right panel shows 6 northern scans, also separated from 
one another by 0'.'5. 



4. Results 

4-.1. Kinematics 

By transforming the spectrum into the velocity space, the 
/ — V map shown in Fig. |^c has been obtained. At a dis- 
tance of 12.2 Mpc, the area covered by the slit of the 
spectrograph projected onto the plane of the galaxy is 
480 pc X 60 pc. The smallest details visible in the I — v- 
map are estimated to be of the order of 0'.'3 which corre- 
sponds to a linear size at the distance of NGC 1386 of only 
~ 18 pc. With this high spatial resolution in our spectra 
it is possible to obtain a detailed picture of the kinematics 
of the NLR in NGC 1386, based on a 2-D echelle spectrum 



The decomposition of the NGC 1386 spectra yielded 
9 individual NLR clouds which show mean velocity dif- 
ferences up to 400kms~^ in [Nii]A6583. The individual 
NLR components, marked by letters, and derived by the 
Gaussian decomposition (see Fig. ^c) are located at dif- 
ferent velocities and spatial locations relative to the nu- 
cleus. The line widths of these clouds vary from about 
90 — 235kms~^ (clouds A, D, F-I) and are accompanied 
by clouds with considerably higher values of FWHM = 
270 - 470kms-i (clouds B, C, and E). The largest and 
most prominent region in this / — u-map is component A 
which is blueshifted with respect to the systemic velocity 
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Hq A6563 








[N II] A6583 






CompoiiBiit 


V [km s 1 


r vv nivi iKiii s 1 


fos. [ J 


r VV xiivi 1 1 


V 1 Km s 1 


r VV niVi [Km b i 


fos. [ J 


r VV nivi 1 1 


at P.A. 16° 


















A 


-125 ± 10 


99 


-1.60 


2.35 


-120 + 10 


89 


-1.60 


2.35 


B 


-200 ± 10 


400 


-0.40 


2.24 


-235 + 10 


471 


-0.35 


2.24 


C 


+200 ± 10 


353 


-0.35 


2.47 


+200 + 10 


377 


-0.33 


2.47 


D 


+75 ± 10 


130 


+1.95 


2.59 


+75 + 10 


153 


+1.95 


2.59 


E 


-145 + 10 


271 


-0.25 


2.59 


-155 + 10 


271 


-0.35 


2.59 


F 


+570 ± 10 


283 


-0.30 


2.35 


+580 + 10 


235 


-0.30 


2.35 


G 


-50 + 10 


118 


-0.85 


2.47 


-35 + 10 


141 


-0.85 


2.47 


H 


-95 + 10 


95 


-0.20 


2.12 


-80 + 10 


118 


-0.20 


2.35 


I 


+35 ± 10 


106 


+0.60 


1.88 


+53 + 10 


113 


+0.50 


1.88 



^ all positional accuracies are of the order of O'.'l. 



by —120 ± lOkms"^ and shows an offset of — 1'.'6 relative 
to the nucleus. 

We have plotted the velocities versus the position in 
Fig. ^c with the derived uncertainties, which show the 
various velocity components for the [N ii] A6583 line. The 
analysis for Ha provided nearly identical results. 

The identified components with their respective veloc- 
ities (relative to the systemic velocity) and their positions 
in our / — ?;-map are listed both for Ha and [N ii] A6583 
in Table |[ The uncertainties arise from the positional ac- 
curacies of our / — w-map. The velocity uncertainties are 
of the order of A w « lOkms^^ and in spatial direction of 

« O'.'l. 

4.2. Morphology of the NLR 

Several emission lines of different ionization levels are cov- 
ered by the wavelength range of our 2-D echelle spectro- 
gram. Therefore, it is possible to search for differences in 
the morphology of the extended line emission as a func- 
tion of the ionization potential. With this method one 
can roughly discriminate NLR components from luminous 
H II regions since H 11 regions are not supposed to be vis- 
ible in the high ionization lines as [Fevii] for instance. 
The mechanism for exciting the gas to temperatures of 
that degree, where highly ionized Fe-lines become appar- 
ent, e.g., [Fevii], can not be achieved with common stellar 
continuum emission as for instance in a massive star burst, 
since starbursts can not supply such high temperatures. 
Therefore, the most likely process is photoionization of 
the gas by a nonstcUar power law continuum provided by 
the central continuum source or by continuum e missio n 
of fast radiative shocks (e.g.,Dopita & Sutherland 1996|) . 
The usual discrimination between Hii regions, AGN-like, 
and starburst galaxy regions is described in the diagnos- 
tic diagrams by plotting the intensity ratios of different 
lOlHl MA cf. Baldwin et al. |l981; 



Luminous H 11 regions in the vicinity of an Seyfert nu- 
cleus are frequently found. These regions are indicative 
of massive starbursts. Sometimes the Hii regions appear 
in a ring structure around the nucleus, which is seen for 
instance in NGC 7469 (DeRobertis & Pogge [1985D and 



Hq 



emission lines (e.g. 

Veilleux & Osterbrock ll987t Allen et al. |l998| ) 



NGC 1068 (Balick & Heckman |1985D . In NGC 1386 there 
are 44 H 11 regions visibl e which were detected by Tsve- 
tanov & Petrosian ( 1995 ) in their survey of Hii regions in 
Seyfert galaxies. These reside in larger radii to the nucleus, 
outside our studied central region. 

The brightest component which we refer to as cloud A, 
was initially suspected to be an ultraluminous H 11 region 
(Rossa 1996), based on the morphology of the 2-D echelle 
line profiles of different excitation levels. The reason that 
NGC 1386 is a low luminosity Seyfert galaxy and a de- 
tection of an apparent ultraluminous H 11 region may lead 
to the conclusion that NGC 1386 is a starburst galaxy. 
However, the detection of the high ionization lines, which 
are characteristic for AGN, rules out that possibility. Fur- 
thermore, the line ratio [Nil] A6583/Ha ~2 rules out the 
possibility even that component A is a H 11 region. It turns 
out that component A is located in the classical region for 
power law excitation in the Seyfert regi me, a ccording to 
diagnostic diagrams (e.g., Baldwin et al. 1981 ). 

The comparison of the morphology of the NLR in the 
Seyfert 2 galaxy NGC 1386 with each of the low ioniza- 
tion lines (e.g.. Ha, [Nil], [Sii]) shows a similar pattern, 
whereby the higher ionization lines like [Fe vii] show a dif- 
ferent morphology. In the low ionization lines the bright 
component A is the dominant source in the I — w-map. 
Comparing this with the velocity field in the [Fe vii] line 
shows that this component is not visible. Here only a cen- 
tral component is appearing which is most likely identi- 
fied with the nucleus of the Seyfert galaxy. Hence, there 
is some evidence that the true nucleus has a much lower 
apparent luminosity than the NLR component A, indicat- 
ing a hidden nucleus. This picture is supported by HST 
imaging studies. It has clearly shown, that nuclear dust 
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Fig. 4. In the left panel the 5 spatial scans (each sep- 
arated by 0'.'65) of the [Fe vii] A 6087 A hne are plot- 
ted and in comparison in the right panel the Ha-|-[Nii] 
complex is shown for the same spatial scans. The ob- 
served wavelength is plotted versus the flux in units of 
10~^^ ergs~^ cm~^ A To better illustrate the differ- 
ences in the Ha-|-[Nii] line scans, we have offset each in- 
dividual scan by 5 X 10"^^ erg s""'^ cm~^ A~^, and in the 
[Fevii] A 6087 A line by 2 x lO'^^ ergs^^ cm-^ A"! 



lanes which are also visible in NGC 1386 are a common 
fea ture w ith typical scale heights of 20-50 pc (Capetti et 
al. 1996 b) which often leads to obscuration of the central 
region of an AGN. Furthermore this is indeed evidenced 
by X-ray measurements, which show that the active nu- 
cleus in NGC 1386 is covered by a Compton-thick (i.e. 
Nh > lO^^cm-^) screen (Maiolino et al. |1998[ ). 

In Fig. ^ we show the spatial variation of profiles of 
the [Nil], Ha lines together with the [Fe vii] A 6087 A line. 



4-3. Identification of individual NLR components 

Next we made a detailed comparison of our reduced 2-D 
echelle spectrogram with WFPC 2 archive images taken 
with the HST. In Fig. ||a we show a superposition of the 
inner part of NGC 1386 of the F658N filter HST image 
with the CASPEC slit. In Fig. |b we show the region (in 
contours) which is covered by the CASPEC slit. In the 
right-hand panel of Fig. ^ the direct comparison of our 
CASPEC spectrum of the [Nil] A 6583 A line as a. I - v~ 
map is shown. The contours are given in logarithmic scale. 
The outer two spectral scans of our echelle spectra are 
severely affected as a result of the flatfielding process as 



described in Sect. 2. We limit our study to the central 
6", where only the outermost regions are to some point 
influenced by the flatfield uncertainties on a small scale. 

Now we can identify the individual NLR components 
seen in our echelle spectrogram with the components in 
the WFPC 2 image. Due to the high dispersion individ- 
ual components of line emission are separated in velocity 
space and can be located spatially with high precision (as 
discussed in Dietrich & Wagner 1998). Therefore, we can 
cross-check the components (or emission knots) in the two 
'images'. This investigation of the morphology can also be 
regarded as an additional method to the direct imaging 
studies in the respective narrow-bands (e.g., [Oiii], Ha). 

In Fig. ^c we have plotted the velocities of the indi- 
vidual NLR components that have been identified in our 
[Nil] A6583A / — ti-map (after Gaussian decomposition) 
as a function of the distance to the nucleus in arcseconds. 
The error bars arise from the uncertainties in the deter- 
mination of the I— and v— coordinates in our / — v map. 

The velocity field of NGC 1386 is not as much dis- 
turbed a s the velocity field in NGC 1068 (Wagner & 
Dietrich 1996) which show s vel ocity differences up to 
1000 kms~^. Weaver et al. ( 1991 ) tried to determine the 
position of the nucleus in their spectroscopic study. Since 
their spectra were taken at medium resolution, the posi- 
tion of the nucleus of NGC 1386 could not be identified 
unambiguously. Most likely a definitive answer would be 
given by polarimetric studies, such as it was done in the 
case of NGC 1068 (Capetti et al. |1996 4 |1997| ), or by as- 
trometry of the radio nucleus. 

The brightest component in our echelle spectrum coin- 
cides with the brightest component in the WFPC 2 image 
of NGC 1386. This bright component A is offset by «95 pc 
with respect to the nucleus. A comparison of the HST 
WFPC 2 image (F658N) with our I - v map of [Nil] A6583 
suggests a correspondence of the components A and G 
with the bright dominating component at ~ — 1'.'5 in the 
HST image of NGC 1386. These two components have a 
similar velocity dispersion of FWHMwlOOkms"^ and are 
blueshifted with respect to the nucleus. Our components 
D and I are associated with the emission knots at « 2'.'0 
and 0'.'5 in the HST image, respectively. The FWHM of D 
and I is comparable to the FWHM of A,G, and H, how- 
ever the components D and I are redshifted with respect to 
the nucleus. The components B,C,E,H might be connected 
with the central bright emission complex at « — 0'.'3 in the 
WFPC 2 image. In contrast to the other components they 
are nearly located along the line of sight and show a ve- 
locity dispersion of the order of FWHM « 350 km s~^ . The 
velocity split of the individual components is of the order 
of 400kms~^ and taking the weak and quite uncertain 
component F into account it rises to Av ~800kms~^. 
Due to the large separation in velocity space from the 
other components we have not plotted component F into 
the I — V map. 
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Fig. 5. Comparison of our CASPEC I — v map with the HST Ha image. The left panel (5 a) shows the central region 
of the WFPC 2 (F658N) image of NGC 1386 (rotated clockwise by 16°) with the CASPEC slit superimposed. The 
second sub-panel (5 b) shows the central region of the F658N image in contours which was covered by the central 6"of 
the slit of our echelle spectrum. The right panel (5c) is the I — v map of the [Nil] A6583A line (shown in logarithmic 
scale) where the velocity relative to the systemic velocity is plotted versus the radius from the center in arcseconds. 
The positions of the individual NLR components (as derived from the Gaussian decomposition) are marked by crosses 
and are accompanied by the respective letters. Due to the large separation in velocity space we have omitted the 
component F in the map. The accurracies of the velocities (10kms~^), and of the positions (O'.'l) are indicated by the 
size of the crosses. 



5. Discussion 

One of the motivations for the investigation of the veloc- 
ity field of the NLR in NGC 1386 is to examine whether 
there are correlations between radio- structures and opti- 
cal emission. In the theoretical picture the deviations from 
a pure rotational component of the velocity are manifested 
in infall or outflow motions. In the case of an outflow 
motion, this often leads to the detection of correspond- 
ing radio structures in the highly resolved VLA maps. 
Those structures are mostly identified with a jet, prop- 
agating through the ambient medium. The interaction of 
the outfiowing radio plasma with the NLR gas will pro- 
duce shocks and this shock gas can be described as an 
ex pandin g cooling co coon around the radio jet (Taylor et 
al. 1992 ; Axon et al. 1998 ). Evidence for this scenario is 



found by Axon et al. ( 1998 ) who mapped spectroscopically 
the inner few arcseconds of NGC 1068 at high resolution. 
They detected a velocity split of the order of 1000 kms~^ 
of individual components which was also found by Diet- 
rich & Wagner ( 1998| ) . The observed velocity difference in 
NGC 1386 can be explained within the framework of the 
model s suggested by Taylor et al. (1992) and Axon et al. 
( 119981) . 

The slit position of NGC 1386 in our inv estiga tion was 
chosen because the study by Mauder et al. (1992) showed 



that the NLR has an elongated structure, with an asym- 
metrical biconical shape. A distinct substructure is visible 
in their speckle images and a few components could be 
resolved in the central part (100 pc). From their speckle 
imaging an opening angle could be derived and the ex- 
tended Ha-emission is mostly visible at a position angle of 
~ 20°. Our identifications of the NLR components are in 
agreement with the observed NLR morphology by Mauder 
et al. ( |199^ ). 

Our simple approach to reconstruct the Ha, [Nil] emis- 
sion line complex using elliptical Gaussian profiles provide 
some evidence for two components of the NLR velocity 
field of NGC1386. Our components can be associated with 
the bright emission line complexes of an HST narrow-band 
image (F658N). The components which might be domi- 
nated by the general galactic rotation can be character- 
ized by FWHMw lOOkms"^ while the components with 
FWHMR:i350kms^^ might be in radial motio n. Th is re- 
sult is in good agreement with Weaver et al. (1991) who 
also provide evidence for a two component velocity field 
of NGC 1386. 

Even the velocity difference of the components with 
larger FWHM is of the order of 400 kms~^ only, there 
might be some interaction of the NLR gas with the out- 
flowing radio plasma. The study by Ulvestad & Wilson 
( 1984 ) revealed a slightly elongated structure in their VLA 
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map (A=6cni) of NGC 1386 at the p.a. 225°. So the radio- 
emission points anti-parallel to the optical emission. How- 
ever, a more recent investigation by Nagar et al. (1999) 
revealed a 300 pc extended feature in p.a. ^ 170° in their 
VLA A3. 6 cm map. 

A further interesting result of our analysis is the de- 
tection of the bright component (cloud A) which is dis- 
placed by — 1'.'6 from the nucleus. In the study by Weaver 
et al. ( 199l[ ) they reported a displacement of the optical 
nucleus by I'.'S from their kinematical studies (kinemat- 
ical center), which is most likely coincident with cloud A 
in our investigation. The true nucleus of NGC 1386 has a 
much less apparent luminosity than cloud A, and hence is 
most likely highly absorbed. This is also evidenced by X- 
ray observations (Maiolino et al. 1998). We conclude that 



component A is excited by photoionization, as is evident 
from the measured line ratios. 

The fact that there is a cone-like geometry of the NLR 
on small scales, hints that absorption plays a role on those 
scales and that effect implies that the ionizing radiation 
can not escape isotropically from the circumnuclear re- 
gions. The observation of cone-like geometries in the cen- 
tral regions of Seyfert galaxies is a common feature. Evans 
et al. ([1991D have shown that the NLR in NGC 1068 is 



cone-like and comprises of several components. The kine- 
matics of the strong components of this cloud ensemble 
have been studied by Dietrich & Wagner ( 199^ ). At least 
~16 components are visible in their I — w-map obtained 
for the [Qui] A5007A-hne down to the scale of 0'.'2 that 
reside in the central 6" region. From their study it was 
inferred that the velocity field of NGC 1068 is driven by 
shocks due to the interaction of the radio jet with the 
ambient medium. 



6. Summary and Conclusions 

We have studied the kinematics and the morphology of 
the NLR in the Seyfert galaxy NGC 1386 with high spatial 
resolution with CASPEC and archival HST images. The 
following results have been obtained. 

— the NLR of NGC 1386 is composed of 9 individual com- 
ponents with 5 components being redshifted and 4 be- 
ing blueshifted with respect to the systemic velocity. 

— the brightest component (cloud A) is blueshifted by 
120±10kms~"'^ and shows an offset of 1'.'6 to the south. 

— this cloud A is identified with the brightest region in 
the HST image, and is excited by photoionization 

— the true nucleus NGC 1386 has a much lower appar- 
ent luminosity, and is possibly hidden by nuclear dust 
lanes. 

— there is evidence for a two component velocity field in 
NGC 1386 

The results show that with this technique of sub- 
arcsecond mapping, which makes use of highly resolved 



2-D spectra, it is possible to identify individual NLR com- 
ponents in the nearest Seyfert galaxies. In this respect it is 
a powerful tool to describe the kinematics and morphology 
of the NLR. 
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